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Abstract 



A conditional symmetry is defined, in the phase-space of a quadratic in velocities 
constrained action, as a simultaneous conformal symmetry of the supermetric and 
the superpotential. It is proven that such a symmetry corresponds to a variational 
(Noether) symmetry. The use of these symmetries as quantum conditions on the 
wave-function entails a kind of selection rule. As an example, the minisuperspace 
model ensuing from a reduction of the Einstein - Hilbert action by considering 
static, spherically symmetric configurations and r as the independent dynamical 
variable, is canonically quantized. The conditional symmetries of this reduced ac- 
tion are used as supplementary conditions on the wave function. Their integrability 
conditions dictate, at a first stage, that only one of the three existing symmetries 
can be consistently imposed. At a second stage one is led to the unique Casimir 
invariant, which is the product of the remaining two, as the only possible second 
condition on ^> . The uniqueness of the dynamical evolution implies the need to 
identify this quadratic integral of motion to the reparametrisation generator. This 
can be achieved by fixing a suitable parametrization of the r-lapse function, ex- 
ploiting the freedom to arbitrarily rescale it. In this particular parametrization the 
measure is chosen to be the determinant of the supermetric. The solutions to the 
combined Wheeler - DeWitt and linear conditional symmetry equations are found 
and seen to depend on the product of the two "scale factors". 

1 Introduction 

In this paper we reexamine the issue of the presence of conditional symmetries 
in minisuperspace constrained systems. The variational symmetries approach for 
classical Bianchi cosmologies has been, to the best of our knowledge, initiated 
in pp. For the case of classical and quantum cosmology, either Bianchi or higher 
derivatives, the first works known to us are [2], [3], [3], [5] and [6] (where conditional 
symmetries are used). Work on the subject has been recently revived, see [TJ, [8], 
[9] and [10]. The method essentially consists in applying the standard theory of 
variational symmetries [12] to the Lagrangian of some minisuperspace model. 
The usual procedure in use is to gauge fix the lapse function to some convenient 
value; then apply the first prolongation in the velocity phase space of a vector in 
the configuration space to this gauge fixed Lagrangian and demand its action to 
be zero. A slightly different approach using the notion of the special projective 
and/or homothetic group is adopted in [T3], [J3] and [To] . 

In our work, we do not gauge fix the lapse function. Our justification is that the 
presence of the lapse keeps intact the full content of the corresponding reparametri- 
sation generator, i.e. the quadratic constraint. This is rather important in various 
ways, not only for the symmetries, as the present work shows, but also in solving 
the classical equations of motion, a task that is greatly facilitated by inserting the 
algebraic solution for N from the quadratic constraint into the other equations of 
motion [16], [IT], [18] and [19]. In the present work we start from a phase space 
point of view and define a conditional symmetry [20] as a simultaneous confor- 
mal field of both the supermetric and the potential. We prove that this definition 
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results in uncovering all variational symmetries. In the Appendix we give an ex- 
plicit example of how the gauge fixing of the lapse function, in conjunction with 
the use of the standard variational condition can lead to a loss of certain existing 
symmetries. 

We give an application of the method to the case of static, spherically symmet- 
ric configurations. A minisuper space Lagrangian for the above family of metrics 
has been given in [21], [22] by considering a 3 + 1 decomposition along the radial 
coordinate r which is taken to be the dynamical variable. In [23], [23] Vakili has 
used the symmetries of the reduced Schwarzschild action in the quantization pro- 
cedure. In the present work we also adopt the quantum analogues of the linear 
integrals of motion as supplementary conditions imposed on the wave function. A 
careful examination of their role and integrability conditions leads us to the unique 
Casimir invariant of their algebra. This invariant is identified to the correspond- 
ing r - reparametrisation generator through a particular redefinition of the lapse 
function. 

The paper is organized as follows: In section 2 we start with some general 
considerations on minisuperspace constrained actions possessing conditional sym- 
metries (see [20 for the case of full pure gravity). Subsequently, in section 3, 
we give the reduced valid Lagrangian reproducing Einstein's field equations for 
static, spherically symmetric configurations, and passing to its Hamiltonian for- 
mulation [25j . we reveal the three conditional symmetries. In section 4 we quantize 
the dynamical system according to Dirac's canonical quantization procedure for 
constrained systems [26] . Finally, some concluding remarks are included in the 
discussion. 

2 General considerations 
2.1 Classical Treatment 

First, let us give some general considerations on conditional symmetries concerning 
reparametrisation invariant actions of finite degrees of freedom. Since we would 
like to address collectively spatially homogeneous and/or point like geometries, we 
adopt as a suitable starting point the following form of the line element: 



ds 2 = ±N(x) 2 dx 2 + g A B(q a (x)) Lof(x l )ujf(x i )dx i dx j i, j = 1,2,3. (2.1) 



The dynamical variable x may either play the role of time in the former or of 
the radial coordinate in the latter case, hence the ± sign in front of dx 2 . Thus, 
N(x), q a (x) are the lapse function and the "scale factor" components, which ought 
to be considered as the dependent dynamical variables of the space-time metric. 
Upon integration of the x l independent variables in the appropriate full action, the 
minisuperspace action A = j Ldx is obtained, with 



L = 



^G ap (q)q a qV-NV(q) 



(2.2) 




a = 1, 



n 
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It is straightforward to verify that the above action retains its form under a 
reparametrisation x = /(£), if one adopts the following transformations of the 
dependent variables q a (x), N(x): 

N(x)^N(x):=N(f(x))f'(x) , q a (x) -> q a {x) := q a (f(x)). (2.3) 

These transformations can be inferred from the underlying geometry (|2.ip . as of 
a kinematical type, dictated by the fact that N(x) appears in the line element 
multiplied by dx while q a {x) multiplied by the differentials of the coordinates 
which are being integrated out of the full action in order to arrive at A. If one 
were presented with only the reduced Lagrangian, then (|2.3j) should be guessed as 
a kind of "gauge" transformation. The above changes of N(x) and q a (x) imply 
that one of these variables can be prescribed at will by suitably choosing the 
freedom encoded in the arbitrariness of /; when one attempts to solve the Euler 
- Lagrange equations, it is both useful and instructive to let the value of N(x) 
be defined by algebraically solving the corresponding constraint equation. If this 
value is substituted into the other equations concerning the q a (x)'s, the rank of the 
system is diminished by one, i.e. one can solve only for n — 1 accelerations. This 
situation is reflected in the presence of two first class constraints in the Hamiltonian 
formulation. 

The momenta corresponding to the configuration variables q a (x), N(x) are 

9L 1 _ n dL . 

vr a := — = — G a0 <f := — r- = 0. (2.4) 

dq a N dN 

According to Dirac's theory 

tt n « (2.5) 

constitutes the primary constraint. The weak equality symbol ~ denotes that the 
corresponding relation can be used only after all the Poisson brackets have been 
calculated. The canonical Hamiltonian is given by 

H:=n 7 f-L = NH C = N G a ^ (q) 7r a tt p + V(q)^j , G^O* = 6&. (2.6) 



The consistency requirement that (|2.5p must be preserved by the x - evolution 
leads to the secondary constraint 

7r N >.= {ir N ,H}n0^H c &0 (2.7) 

and the algorithm is terminated since {T~L C , H} = which is stronger than what 
is actually needed, i.e. ~ 0. The two constraints are first class on account of the 
Poisson bracket {ttn^c} = 0. They represent the x - reparametrisation invariance 
of the action and reveal N(x) as a Lagrange multiplier (as pointed out above) not 
defined by the dynamical evolution equations involving the accelerations of the 
n — 1 g Q (x)'s. 

On the part of the configuration space spanned by q a (x), one can define a con- 
ditional symmetry generated by a vector field £ which is a simultaneous conformal 
field of the metric G a/3 (q) and the potential V(q). This leads to the following: 
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Definition. A vector field £ Q on the configuration space, generates a conditional 
symmetry if 

C^ = ^(q)G^ , C 6 V(q) = </>(q)V(q) (2.8) 

where stands for the Lie derivative operator acting on the corresponding 
geometrical objects: 

:=Gfe-G^^ p -G a P^ p , C ( V(q):=V(q) >p e- (2.9) 

To each of the existing conditional symmetries one can correspond a phase 
space quantity 

Qj := C?7T a . (2.10) 

As one can easily check, the Poisson brackets of (|2.10p with the Hamiltonian vanish 
on the constrained surface on account of (|2,7|) : 

{Qi,NH c } = -N \^(c^G a ^-K a T:p+ C^v\ = -N fafa) U c « 0. (2.11) 

Therefore, by virtue of (|2.11|) the quantities (|2.1U|) are constants of motion 

Qi = ki. (2.12) 

As it is well known, the totality of these charges forms a Lie algebra characterized 
by some structure constants, say Cfj 

{Ql,Qj} = CfjQ M . (2.13) 

The equation {Qi, H} ~ is the definition for conditional symmetries used in 
[20] . Our definition of conditional symmetries (|2.8p is of course equivalent to the 
aforesaid definition, so it is a matter of taste which one should be preferred over 
the other. Moreover, (|2.8p are somewhat different than the conditions used in the 
earlier works cited in the introduction, in the sense that we allow the £/'s to be 
conformal fields of the supermetric and not Killing fields. This is justified by the 
weak vanishing of the superhamiltonian H resulting from (|2.7p . 

A further important property emanating from this definition is the fact that 
under a scaling of G°P and V{x) by, say ui(q), the conserved charges (|2.10p retain 
their nature and, more usefully, remain form invariant. To see this, let us act with 
the same on G al3 = u(q) G afi and V = u(q) V(q): 

£ e (Lo(q)V) = u:{q) C^V + C^(q) V = U + u{q)V =: <%) (u(q) V) 

(2-14) 

and likewise for G a P. The above equation shows that if we form the new quantities 
Qj = £f it a , these remain constants of motion, since their Poisson bracket with the 
He still vanishes weakly. 

The scaling by u> can be achieved by rescaling N, i.e. iV = u(q)N, as it can 
be seen by the form of (|2.6p . A particularly interesting and useful choice of uj(q) 
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is the one that makes $(q) equal to zero. This uj(q) can be evaluated with the aid 
of the second of (ESI): 



c^v = <j>{q)v ^ev a = m v => m = (2.15) 

Use of this 4> into the definition of <&(q) leads to a linear partial differential equation 
for oj which can always be solved. 

An obvious solution is u) = XV" 1 (with A being a constant) which makes the po- 
tential V constant. Furthermore, in this case the £/'s are turned into Killing fields 
of G al3 rather than conformal fields. It is in this particular parametrization for the 
lapse that a connection can be established with the usual variational (Noether) 
symmetries approach: The Killing fields £/, with the appropriate prolongation, 
provide the fields Xj which are the generators of the corresponding symmetries, 
i.e. they satisfy pr^Xj(L) + LDivS/ = DivFj [12] (with Sj being the coefficient 
of the partial derivative of the independent variable x and Fj(q a , x) being an arbi- 
trary function). This is rather important, since if one applies the above condition 
to another lapse parametrization in which the potential is not constant and the £/'s 
are not Killing vector fields, one is bound to loose some symmetries, if not all as 
the Appendix shows for the case of the specific example of section 3. On the other 
hand, our definition of symmetries (|2.8p . based on the phase space and inspired 
by the weak vanishing of the superhamiltonian, can be applied to any conformal 
gauge. Therefore we have the following theorem: 

Theorem. For the case of constrained systems described by f)2.2[) , the con- 
ditional symmetries are equivalent to Noether symmetries. 

If one wished, one could also work in the velocity phase space and in an arbitrary 
parametrization of N(x). The price one would have to pay is to accordingly modify 
the variational symmetry condition by allowing the right hand side to be expressed 
as a multiple of the constraint Euler - Lagrange equation with respect to N. Let 
us start with 

X:=^(q)^-, pr%:=I + ^MA (2.16) 
"» KH> dq a 1 dx dq a y 1 

Definition. A vector field X = £, a (q)d a generates a conditional symmetry if 

pr^X{L) = (j){q) N EL(N) (2.17) 
where EL(N) is the Euler- Lagrange equation with respect to N . 

Using the above definition we can easily compute the action of pr^X on (|2.2p 

pr (1) X(L) = <p{q) N EL(N) 
± (CtGap) q a q? - NC^V = -<P(q) N (-L + V \ ^ 

± (CtGap + 4>{q) G a p)q a q P ~N (C^V- </>(q) V) = V<f (2.18) 
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which, in order to be an identity for all q^ leads to the conditions (j2.8|) thus 
proving the equivalence of the two definitions. This can also be visualized as an 
enlargement of the relevant vector X by adding the term —cp(q)N-^, making 
manifest the freedom to arbitrary rescale N. 



2.2 Quantum Treatment 

For the canonical quantization of the above system, we adopt the point of view 
of promoting to operators not only 7i c , but also the generators of the conditional 
symmetries Qj. In the Schrodinger picture the first class constraints become the 
quantum conditions: 

d 

TC N V(q, N) := -i— *(g, N) = * = 



H r ^ = 



V = (2.19) 



where u is a suitable measure which must transform as a scalar density; under 
u the quantum quadratic constraint (|2.19p is Hermitian. This form is the most 
general second order scalar operator, except for a possible multiple of the Ricci 
scalar, 1Z, corresponding to G a p. 

A natural scalar density emerging from the geometric structure of the configu- 



ration space is det G a p\ = y/\G\ . It is reasonable to limit the possible measures 
to be inferred from the classical geometrical structure, as n = f(y/\G\). How- 
ever, this choice of \x in conjunction to the scalar density transformation law for 
\x and \J\G\ leads to the equation p,(y/\G\J) = fi(y/\G\)J, which must hold for 



any transformation q a = q a (q^) with J = det The only solution is therefore 

H = e^ q ' y/\G\, with A any scalar. Thus, one must essentially select the Laplacian 
based on the configuration space metric, if X(q a ) = const., and the addition should 
then be exactly 4 ^"~ 2 ^ 1Z in order for the operator to be conformally covariant |27| . 

The conditional symmetries (|2.10p . f|2. 12[) must also be turned into Hermitian 
operators (see [28J), as : 

Qi*:=-^(ut;?d a + d a utf)y (2.20a) 
Qjty = Kl y. (2.20b) 

If we bear in mind that the classical quantities Qi generate symmetries, it is 
mandatory, in order to retain their geometrical character, to demand that they act 
as derivatives on the wave function ^ . In other words, we must have the momentum 
operators acting on the far right side, which translates into the condition -^— t — — = 
0. This equation can be satisfied in the particular conformal parametrization 
N = y, since then 

« . W A >i^ff) . w + i%a . „ ^ w . „. ,, 21) 

\i exp (A) \G\ l l l 
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The term £f. a is zero if evaluated with the metric G a p in which the £/'s are Killing 
fields. 

The maximum number of Qi's is ■ If the actual existing symmetries are 

more than or equal to n, the above equation fixes A to a constant value. For a 
lesser number of symmetries, some non trivial X(q) may be allowed. 

Upon quantization it is reasonable to assume that a quantum algebra isomor- 
phic to the classical one holds: 



which turns (|2. 13j) into 



<■■■>- x 



(2.22) 



Qi, Q, 



±hC M u Q 



IJ WM 



Qi,Q^ 



±C M U Q 



I J HM 



(2.23) 



where in the last equality we used units in which h = 1. At first sight equation 
(|2.23|) in conjunction with ()2.20p seems to impose certain restrictions on the form 
of the £'s, the constants ki and an initially arbitrary measure /x(g a ) (although we 
have already committed ourselves to fi oc ^/|G|). In order to begin exploring these 
restrictions, we act with f|2.23|) on ^(q a ) and we use (|2.20b[) 



Qi,Qj 



y(q a ) = ±C M jjQ M *(q 



(kj kj - kj « z ) *(g Q ) = i C M U K M *(q a ) 



C M IJ KM = 



(2.24) 



The above equation must be interpreted in a recursive and exhaustive manner: It 
must be checked for all subalgebras since they will correspond to different admis- 
sible choices of subsets of the Q/'s, / C {1, . . . , |, The number of essential 



constants for the underlying geometry ()2.ip can provide us with a lower bound of 
the dimension of the subalgebras to be selected: If (|2.ip depends on some constants 
K\, K2, one can always try to find a coordinate transformation of x, x l that ab- 
sorbs some (or all) of the constants. If a constant cannot be absorbed by such a 
transformation we call this constant essential. From this definition it is apparent 
that the constants not appearing in (|2.24p are eligible candidates for represent- 
ing the geometry. Such constants will belong to the centre of the corresponding 
subalgebra. 

Further restrictions, induced by (|2.23|) . could come to life when we expand this 
equation taking into account (|2.20ap . The left hand side of (|2.23p . then reads 



l.h.s 



2/x(g* 



{»{q v )Hi a ) a -Hfd p + 



+ 



2^) \»Wtj Jtfl 



tj p fy, 



(2.25) 



The first commutator of (|2.25p can be evaluated quite easily if we expand (|2.13p 



e/ ^ - e/ tj a P = c m jj e M Q , (2.26) 

while the last two commutators of (|2.25p combine to 

2^(^) ^ j/3 ' Q ~~ ^ jCl ^ //3 ' Q ) + 2 ^ j/3 ' Q / ? ~~ ^"Z 3 ) ' (2.27) 
The right hand side of (|2.23p reads 

r.h.s. = ±C M U (-±i M a d a + KM°), a 

= ^ 7 5 q + ( M (<f) £ m ") q . (2.28) 

Gathering the above results we get the integrability condition 

+C M /;(M«^M°), a = (2.29) 

and using once more (|2.26p we arrive at 

C M U V£ M + [Ci VO - Cj V£/) = 0, (2.30) 

where V£/ = d a ^f. Surprisingly enough, the above relation is nothing but an 
identity that can be derived by taking the divergence of ()2.26p . leaving us with 
only (gjUp . 

Finally, let us point out that if someone wished to investigate the possibility 
that the quantum algebra is not isomorphic to the classical algebra, thus relaxing 
(|2.23p and just evaluating the commutator on the wave function <J>(g a ), then one 
would obtain a combined relation of the form 

i C M jjK M - \ (A^VO - £ V& + C M jj V£ A f) = (2.31) 

which, however, is again reduced to (|2.24p . This establishes (|2.24p as the only 
constraint on the classical quantities, even without the assumption of the quantum 
algebra being isomorphic to the classical algebra. 

These constraints among the k^s prohibit the simultaneous realization of all 
Qls as eigen-operators. Usually one adopts the maximal Abelian subgroup along 
with the Casimir invariant as consistent conditions on the wave function. However, 
choices of non- Abelian subgroups may also be consistent [29]. Further investigation 
can be made for particular systems as in the example below. 



9 



3 Hamiltonian formulation of static, spheri- 
cally symmetric geometries 

Our starting point is the static, spherically symmetric line element 

ds 2 = -a{rfdt 2 + n{rfdr 2 + b{rf{d6 2 + sin 2 dcp 2 ), (3.1) 

where we have changed the notation of the lapse function from N(x) to n(r). 

In the usual ADM 3+1 decomposition one foliates space-time in t-hypersurfaces 
and the coefficient of dt 2 is the lapse function. Here we adopt a 3 + 1 foliation in 
the r coordinate and therefore the role of the lapse is attributed to the coefficient 
of dr 2 in the line element. Thus we consider n(r) in (|3.1|) to be the r-lapse function 
while a(r) and b(r) are the dependent "dynamical" variables on the r-hypersurface. 
Previous authors use also a "shift" term 2B(r) dr dt. However, this is not relevant 
for the following discussion, as B(r) does not enter the Einstein tensor while it can 
be absorbed by a time redefinition of the form t = t + J dr. 

The Einstein - Hilbert action Ae~h = J \/—gRd 4 x for the geometries (|3.ip 
leads to the reduced action A = j Ldr with the following Lagrange function 
L(a, b, a' , b',n): 

Aba'b' 2ab' 2 , . 

L = 2an + + , (3.2 

n n 

where ' denotes differentiation with respect to the spatial coordinate r. It is easy 
to verify that the Euler - Lagrange equations obtained from (|3.2p are identical to 
Einstein's equations Q^ u = for the line element (|3.ip . 

In order to proceed with the Hamiltonian formalism we calculate the conjugate 
momenta, 



(3.3) 





dC 






TTn = 


dn' 


= 0, 






dC 


Abb 1 




7T a = 


da 1 


3 

n 






dC 


Aba' 


Aab' 


7T b = 


db' 


n 


h 

n 



Obviously, 7r ra is a primary constraint. The Legendre transformation leads to the 
Hamiltonian 

H = nn c , 

where 

V - 9n a ^ a + 7Ta7Tb (Id) 
^--2a-^ + — . (3.4) 

The preservation of the primary constraint ir n in the r— evolution, i.e. 

< = K,#}«o, 

leads to the secondary constraint 

H c « 0. (3.5) 
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The minisuperspace metric inferred from (|3,4p is 

G Q/3 = . (3.6) 

V & o) 

Our definition (|2.8p is fulfilled by the following three conformal Killing fields of 
both G a p and the potential V = —2a: 

(1 = (-„,„>, 6 =(i,0), 6=(-i,l) (3-7) 
which, contracted with (7r a ,7T&), provide us with the three integrals of motion: 

Q x = -a7T a + bTT b , Q 2 = -^r, Q 3 = — ^+7T6- (3.8) 

a o 2 b 

We calculate the Poisson brackets of these conserved quantities with the canon- 
ical Hamiltonian H and the Poisson algebra that they satisfy: 

Tl Tl 

{Q 1 ,H} = nH c , {Q2,H} = -—H C , {Qs,H} = —n c , (3.9a) 

a z o lb 

{Qi,Q 3 } = Q 3 , {Q2,Qi} = Q 2 , {Qs,Q2} = o. (3.9b) 

As expected from the discussion of the general case (|2.1ip . the Poisson brackets 
(|3.9ap are weakly vanishing on the constraint surface % c ~ and therefore the 
three Qi's are constants of motion. 

At this point it is interesting, and useful for what follows in the quantization, 
to adopt a new parametrization of the lapse n(r) = 2a(r) which makes the po- 
tential constant as explained in the previous section. The Lagrangian and the 
corresponding Hamiltonian are now given by 

- Saba'b' Aa 2 b' 2 - - 1 

L = n+ = + — = — , H = nn c = n—n c . (3.10 

n n 2a 

If the value of n specified by the constrained equation is substituted into the Euler 
- Lagrange equations for a(r) and b(r) the system can be solved for only one 
acceleration, say a"(r), and the general solution of the entire system is: 



n(r)=2cb'(r) , a (r) = c^l-^, (3.11) 

where the constants of integration have been rearranged so that the ensuing line 
clement 

ds 2 = -c 2 (l- ^r) dt 2 +( 1 - ^r) 1 b'(r) 2 dr 2 +b(r) 2 (de 2 + S m 2 9d^ 2 ) (3.12) 

V o{r)J V KnJ 

bears the closest possible resemblance to the standard form of the Schwarzschild 
metric, while the presence of the arbitrary function b{r) reflects the r reparametri- 
sation covariance of the system. 
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The new supermetric is now given by 



'W Sab\ 

G ap = I , (3.13) 

8^6 / 

and the fields (|3.7|) are now turned into Killing fields of G a ^; they are also, trivially, 
symmetries of the constant potential V = 1 as well. The algebra satisfied by the 
new quantities Qi = ^fn a and % c can be easily seen to be 

{Q 1 ,n c } = 0, {Q 2 ,H c } = 0, {Q 3 ,H c } = 0. (3.14) 

Therefore the quantities Qj still remain constants of motion. It is in this conformal 
gauge that one can explicitly verify that the prolonged £/'s, Xj, are satisfying a 
form of the standard condition of variational symmetries pr^Xj(L) = 0, thus 
ensuring that they are generators of Noether symmetries. It is important to note 
that if one had inserted the Xj's into the full relation pr^'X(L) + LDivH = DrvF 
with L given by (|3.2p the result would be negative. This situation is clearly depicted 
in the Appendix. 

If we invert the definitions for Tt a and use (|3.1ip in Qi we calculate the integrals' 
values on the solution space: 

Q l =Q 1 =AcM , Q 2 = Q 2 = - , Q 3 = Q 3 = 2c (3.15) 

c 

As it is known the Schwarzschild solution involves only one essential constant, 
the mass M. The second constant appearing in the solution (|3.1ip . (|3.12|) can be 
seen to be absorbable by a scaling of the time coordinate t — > -, allowing to set 
c = 1 but not c = 0. So, on the solution space we can set Q\ = AM, Q2 = 4, 
Q3 = 2. It is noteworthy that the values of the last two integrals of motion can 
be changed at will and that they do not involve the essential parameter M which 
characterizes the geometry. The above argument relies on the underlying geometry. 

But what if we were deprived of the line element and we were just given the 
dynamical system (|3.10p ? How could we differentiate between the constants K\ 
and «2j K 3? Interestingly enough, there is an argument that leads to a distinction 
between them: The crucial observation is that Qi h as a vanishing Poisson bracket 
with each kinetic term of H, while Q2, Q3 need the entire kinetic term in order 
to produce vanishing Poisson brackets. This fact is reflected in the following prop- 
erty concerning the one-parameter family of canonical transformations generated 
by the charge Q\ (see [30], [31] for the generalization of Noether symmetries for 
constrained systems and for Noether's theorem in phase space), 

a e x a, Tt a — > e~ X 7t a , b — > e~ x b, itf, — > e X 7tb (3.16) 

Under such transformations Qi remains, of course, unchanged, while Q2, Q3 are 
scaled by e~ x and e x respectively. One can thus use the freedom of A to arbitrarily 
change the values K3, but not k±. Furthermore, the Hamiltonian T-L c remains, 
due to the particular scaling of Q2, Q3, unchanged. 
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If we return to the phase space, we can write 

023 = Q2Q3 = -TTTo^a + —T^ct^fi (3.17) 

2 cr ab 

with the relevant Poisson bracket algebra now becoming 

{Q 23 ,^ c } = (3.18a) 
{Q/,0 23 } = 0, 7=1,2,3. (3.18b) 

As expected Q23 is a quadratic integral of motion. From a group theoretical view, 
Q23 is an element of the centre of the universal enveloping algebra (uea) generated 
by Q/s, i.e. it is the Casimir invariant. The Hamiltonian T-L c belongs also to the 
centre of uea (it commutes with all Q/'s), thus it can only differ from Q23 by an 
additive and/or a multiplicative constant. It is an easy matter for one to check 
that indeed 

Q 23 = 8(^ c + 1). (3.19) 

To sum up, we have constructed a gauge independent, quadratic in the mo- 
menta, integral of motion which commutes with the only integral of motion that 
entangles the sole essential constant of the Schwarzschild solution. In the next 
section, and in order to proceed with the quantization, we will rely on these two 
quantities. 



4 Quantization 

In order to quantize our system, we must turn into operators H c = g (Q23 — 8) , 
Qi, Q2 and Q3 (hereafter, for the shake of simplicity, we will omit the bars from the 
symbols of the corresponding operators). The corresponding quantum operators 
can be inferred from (|2.19|) . (|2.2Ua|) . (|3.13|) and (|3.7p . On account of our previous 
general discussion, the measure entering the quantum operators ought to be taken 

as n(a,b) = Xy/ det\G a p\ at ab. The constancy of A is forced by the combined 
requirement that the Q/'s must be realized as Hermitian operators and at the same 
time retain their classical geometrical character by acting as derivatives. Thus the 
extra term ifd a \ must vanish for all 7 = 1,2, 3, which leads to a constant A. Two 
further arguments in favor of this choice of measure are: 

• The fact that the quantum analogue of the algebra (|3.18bp is made isomorphic 
to the classical, i.e. 

[Qi,Q2 3 ]F(a,b) = 0, 7=1,2,3 for any F(a, b) (4.1) 

a fact that is highly non trivial, since it depends on the choice of both the 
factor ordering and the measure. 

• At the classical level, the only linear integral of motion involving the es- 
sential constant is Q\. If we seek the functions on the configuration space 
which are invariant under the point transformations generated by Q\ we find 
{Q 1 ,f( a ,b)} = 0^f(a,b) = f(ab). 
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The above arguments lead to the following linear operators corresponding to the 
elements of the classical algebra, the Casimir invariant and the Hamiltonian: 



Qi = -±(bd b -ad a ) 



(4.2) 




(4.3) 



Q 8 = -i(ft-^0a) 



(4.4) 



^23 = T5 d a d a 7 9 a d b + 



1 



(4.5) 



2a b 2 



n c = \ (q 23 - s) 



(4.6) 



It is an easy task to check that these operators satisfy not only the relations 



for any test function F(a, b), but also (|4,ip as well. 

Due to the constraint condition (|2.24p . applied to the specific structure con- 
stants inferred from (|3.9bp . we conclude that only the eigenvalue k,\ is free and K2, 

must necessarily be zero. However, the latter is impossible since on the classical 
solution space f)3. 15j) hold. The reasoning given below (|2.24[) results in the need to 
consider the two- and/or one- dimensional subalgebras. The investigation of these 
cases can be easily carried out. 

As far as the 2d subalgebras are concerned the results obtained are briefly the 
following: 

(a) For the two non Abelian subgroups either K2 or K3 are forced to be zero, 
something that is inconsistent with their classical values. 

(b) For the Abelian subgroup, the two linear equations lead to the solutions 
ty(a, b) = A exp (§(ft2 a 2 b + 2ks 6)), where A is constant and the quadratic 
constraint enforces the restriction K2 K3 = 8. It is of course doubtful if one 
can accept such a wave function to represent the geometry, knowing that it 
does not contain the essential constant M. However, one could interpret it 
as plane waves representing the limiting flat space-time M = 0. 

The one-dimensional subalgebras spanned by Q2, Q3 give solutions which are spe- 
cial cases of the solution described in b), as expected since they commute. 

Consequently, the only possibility is to adopt H c and Q\ as conditions on the 
wave function. This is indeed possible, since they commute with each other and 
therefore can be considered as physical quantities on the phase space that can 
be measured "simultaneously" (our dynamical parameter is the distance r). The 
ensuing eigenvalue equations are 



[Qi,Qj]F{a,b) = ±C K u Q K F(a,b) 



(4.7) 



= => a d a d a ^-26 d a d b * + d a * - 16 a b 2 = 
Qi* = «i * =>• i (-ad a $> + bd b ^>) = 



(4.8b) 
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The solution of the linear partial differential equation (|4.8bp is 



$(o,6) = a iKl S(ab). (4.9) 

If we insert the above solution into the Hamiltonian constraint (|4.8ap we arrive at 
the following ordinary differential equation for S(u) (u = ab): 

u 2 S"(u) + uS'(u) + (k? + 16 u 2 ) S{u) = (4.10) 

which has the general solution 

S(u) = ci J iKl {Au) + c 2 Y iKl (4u), (4.11) 

in terms of the Bessel functions of imaginary order. 

In order to gain some insight on the normalizability of the formal probability, 
instead of these Bessel functions and because of their imaginary order, we can use 
the functions F± K1 (4u) and G± Kl (4u) defined in [32] through the Hankel functions 
h£\u) = J^{u) + lY^u) and H^{u) = J^u) - tY^u), \i G C. Thus, the 
solution can be written as 

S(u) = c 1 F iK1 (4u) + c 2 G iKl (Au) (4.12) 

with 

F iK1 (4u) = \ [e-^l 2 fl« (4«) + e^' 2 flg> (4«)) (4.13) 

G iKl (4n) = 1 (e~™ /2 ^£(4n) - e^' 2 H^u)) . (4.14) 

These functions are linearly independent solutions of (|4.10p and have the following 
properties: a) when u G (0, +oo) they are real, b) they are oscillatory with a phase 
difference of \ and c) when both u, and/or k\ tend to zero, F± R1 (4u) tends to 1, 
while G± K1 (4u) becomes infinite. 

The final form of the wave function ^(a,b) is 

*(o,6) = a iKl S(ab), (4.15) 

so we can define a probability density of the form 

n(ab)V*(a,b)V(a,b) oc u S*(u) S(u). (4.16) 

5 Discussion 

In this paper we give a definition of conditional symmetries (|2.8|) in terms of 
simultaneous conformal Killing fields £/ of the minisuperspace metric G a/3 and 
the potential V. We prove that these symmetries remain form invariant under a 
rescaling of the lapse function N = uj(q)N. Accordingly, we observe that there is a 
special scaling of the lapse by the inverse of the potential u = y , in which the form 
invariant £/'s become Killing vector fields of G a ^ . In this conformal gauge we show 
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that these £/'s are indeed equivalent to the variational (Noether) symmetries one 
would have found if one had used the standard definition of variational symmetries 
on the specially scaled Lagrangian L = ^G a/ 3q a q 13 — N. We also show that it 
is possible to find the same variational (Noether) symmetries in the arbitrary 
parametrization L = ™ G a ^q a cf — N V if (|2,17p is used in the configuration 
space. 

As a result, we propose a method to quantize the minisuperspace actions which 
are described by the singular Lagrangian (|2.2p . The steps of this procedure are: 

1. Go over to the Hamiltonian H. 

2. Calculate the Noether symmetries as Killing fields of the metric G a ^{q) and, 
trivially, symmetries of the constant potential. 

3. Identify the essential constants of the metric via the arguments following 
(|2.24p . Promote the allowed Q/s to operators according to (|2.20p with \i = 
|G|V2. 

4. Promote the Hamiltonian constrain T-L c (which is a linear function of the 
Casimir invariant) via (|2.19p to a Hermitian operator acting on the wave 
function. 

In sections 3 and 4 we present an application of the above method for the case 
of static, spherically symmetric geometries. First, we begin from the Lagrangian 
(|3.2p emanating from the line element (|3.ip . We find the simultaneous conformal 
Killing fields (|3.7p of the supermetric and the potential, which define the three 
conserved charges (|3.8p . The unique Casimir invariant of their algebra is (J23- In 
order to make it numerically proportional to the kinetic part of the Hamiltonian 
we are led to (|3.1U|) . 

• In [3J Capozziello and Lambiase use the standard Noether symmetries ap- 
proach for the regular system obtained by gauge fixing the lapse function 
occurring in a singular minisuperspace Lagrangian. Furthermore, they pro- 
pose the use as quantum operators of as many of the symmetries as they can 
be simultaneously brought into normal form by a single coordinate transfor- 
mation of the configuration space variables, so that they become manifestly 
cyclic. They thus, effectively, invoke the maximal Abelian subgroup as the 
relevant tool for quantization. Our perspective is quite different, since we 
start with a singular Lagrangian (|2.2p and we use only those symmetries 
that are allowed by the condition (|2.24p . In this sense, we have refined their 
method and their search for "... a criterion by which the Hartle point of view 
can be recovered without arbitrariness" ; the selection rule (|2.24p is exactly of 
that nature. An other worth-emphasizing point is that the use of the Noether 
symmetries must be restricted by the requirement that they correspond to 
essential constants of the underlying geometry; in our example, the use of 
the maximal Abelian subgroup leads to the marginally acceptable plane- wave 
solutions not containing the classical geometry's essential constant M. 

• Our specific example (described in sections 3 and 4) has also been the sub- 
ject of [23]. In that work, Vakili finds two of the Q/'s but then he uses a 
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linear combination of them in order to reproduce the essential constant M 
of the Schwarzchild metric. He thus reaches to the unique acceptable linear 
quantum operator equivalent to our Q\. The clever choice of the lapse func- 
tion, his equation (4), along with the somewhat unorthodox choice of factor 
ordering for the operators (see below his equation (55)) leads essentially to 
a constant potential (his equation (53)) and to the Laplacian operator (his 
equation (56)). As a result, the solution spaces found both by us and Vakili 
essentially coincide. Of course, our general theory constitutes a systematic 
explanation of the various choices of his work. 

• In [33J Jizba and Pons use Noether symmetries in order to transform a regular 
Lagrangian into a singular one: they promote, at the classical level, the con- 
stants of motion to constraints by adding them, with appropriate Lagrange 
multipliers, to the Lagrangian. The most natural step after this, in order to 
quantize the theory dla Dirac, is to promote the constants of motion to oper- 
ators acting on the wave function. This is exactly our way of thinking since 
we already start with a singular Lagrangian and we apply the appropriate 
Qi on the wave function. 

A further point that we would like to stress is the somewhat unexpected result 
that the quantum algebra (|2.24p of the linear operators, turns to be isomorphic 
to the classical one (|2.13|) by virtue of the form of the operators (|2.20p and the 
relation among the £'s (|2.26[) inferred from (|2.13p . 

We plan to return with an exhaustive list of applications to all 2d and/or 3d 
configuration minisuperspaces, emanating from the appropriate spacetime geome- 
tries, i.e. Bianchi types. 
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A Appendix: Calculation of variational sym- 
metries 

As already stated in the main text, if one searched for the generators Xj of varia- 
tional symmetries of (|3.2|) by using the relation 

pr^Xj(L) = 0, (A.l) 

none of the three £'s would emerge. On the other hand, condition (jA.ip works for 
the reparametrised Lagrangian (|3.10|) and reveals the three conditional symmetries. 
This is due to the fact that we are dealing with a singular Lagrangian. In the first 
case, the Qi = (,fir a are integrals of motion on the reduced phase space H c « 0, 
while in the latter, {Qi, H} = 0. 

In order to calculate the right X/'s in any reparametrisation we use the def- 
inition (|2.17p in which there exists an extra term proportional to the quadratic 
constraint equation involving the velocities 

pr^X^L) = <f>(a,b)nEL(n) (A.2) 

In the case of the Schwarzschild example we have 

X :=rj a (a,b)—+r] b (a,b)— (A.3) 

r (1) X •= X + — — + — — (A 4) 

dr da' dr da' 

and by substituting (|3.2p inside (|A.2p it is easy to calculate that 

d , d v Id ad d 

X i = - a d-a +b db> * 2 = ^' X * = -2bd-a + 0b> (A - 5) 

which are exactly the £j's of ([3.7]) . The corresponding functions 4>i that multiply 
the Euler - Lagrange equation EL(n) = ^ are 



1, 2 = ^, ^3 = -^ (A.0, 



and are equal to the multiplying factors in (|3.9ap divided by —n. 

As a concluding remark it is reasonable to say that, it is of utmost importance 
not to fix the gauge, i.e. set n = 1, since the presence of the quadratic constraint 
is needed in order to acquire all the variational symmetries regarding a singular 
Lagrangian. This, of course, does not mean that one is prohibited to select a 
gauge for n; one must simply remember to take into account the gauge fixed form 
of the constraint equation as in ()A.2p . The extra term in (|A.2p can be interpreted 
as a component of a more general generator, Xj = —cfiind n + Xj, expressing 
not only transformations of the reduced configuration space (a, b) as X/'s do, but 
transformations over the full space (n, a, b). This does not add any term of the 
form tJ^- to the prolonged vector, since the Lagrangian is free of n'. 
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